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ABSTRACT 

A simple algorithm is presented for the calculation of the projected line-of-sight ve- 
locity profiles (VPs) of non-rotating anisotropic spherical dynamical models with a 
phase-space distribution function of the Osipkov-Merritt type. The velocity distribu- 
tion in these models is isotropic inside the anisotropy radius r a and becomes increas- 
ingly radially biased at larger radii. VP shape parameters are presented for a family 
of models in which the luminous mass density has a power-law cusp pi, oc 7 1-7 at small 
radii and a power-law fall-off pi, oc r~ 4 at large radii. Self-consistent models and 
models in which the luminous matter is embedded in a dark halo are discussed. The 
effects of changes in the cusp slope 7 and in the anisotropy radius r a are documented, 
and the area in the (7, T , a )-plane that contains physical models is delineated. The 
shapes of the VPs of the models show a considerable (and observable) variation with 
projected galactocentric radius. These models will be useful for interpreting the data 
on the VP shapes of elliptical galaxies that are now becoming available. 

Key words: galaxies: elliptical - galaxies: kinematics and dynamics - galaxies: struc- 
ture - line: profiles. 



1 INTRODUCTION 

The stellar kinematics of elliptical galaxies provides informa- 
tion on the process of galaxy formation and on the presence 
and properties of dark matter components, such as a dark 
halo or a central black hole. Traditionally, only the mean 
stellar line-of-sight velocity and velocity dispersion could be 
determined from observed galaxy spectra. Only recently 
(e.g., Franx & Illingworth 1988) has it become possible to 
derive from observations the full projected line-of-sight ve- 
locity distribution, henceforth referred to as the velocity pro- 
file (VP). A variety of data analysis techniques is now avail- 
able for this purpose (Bender 1990; Rix & White 1992; van 
der Marel & Franx 1993; Kuijken & Merrifield 1993; Saha 
& Williams 1994). These have been used to determine VP 
shapes for a significant number of galaxies (e.g., Bender, 
Saglia & Gerhard 1994; van der Marel et al. 1994a; Car- 
ollo et al. 1995). Stellar dynamical models are required to 
interpret such observations. 

Combined photometric and kinematic data generally do 
not uniquely constrain the stellar phase-space distribution 
function (DF) / (e.g., Dejonghe 1987). The conventional 
approach is therefore to calculate observable properties for 
special classes of DF's, and compare these to data. In this 
paper we restrict ourselves to non-rotating spherical sys- 
tems. Therefore the DF depends on the energy and squared 
angular momentum per unit mass, / = f(E, £ 2 ) (Binney & 
Tremaine 1987, hereafter BT). Many DFs of this type have 



been discussed in the literature, some chosen for mathemati- 
cal convenience, others for physical plausibility (e.g., Michie 
1963; Osipkov 1979, hereafter 079; Merritt 1985a,b, here- 
after M85; Stiavelli & Bertin 1985; Cuddeford 1991; Gerhard 
1991; Louis 1993). 

One constraint to be satisfied by any realistic DF is that 
it matches the observed surface brightness profiles of ellipti- 
cal galaxies. De Vaucouleurs' (1948) R 1 ^ 1 law generally pro- 
vides a remarkably good fit to ground-based observations. 
However, the corresponding mass density and gravitational 
potential cannot be evaluated analytically, which compli- 
cates the construction of dynamical models. For this rea- 
son, various simple analytical potential-density pairs have 
been proposed that in projection approximate the R 1 ^ 1 law. 
Many of the popular models have a constant mass density 
core. Jaffe (1983) and Hernquist (1990, hereafter H90) re- 
alized that such models do not fit the R 1 ^ 1 law at small 
radii, and proposed simple analytical potential-density pairs 
with a central mass density cusp. In their models the mass 
density inside the half-mass radius varies as p oc r~ 2 and 
p oc respectively. Recently, observations with the Hub- 
ble Space Telescope have revealed that essentially all early- 
type galaxies have central surface brightness cusps, with a 
large variety of cusp slopes (e.g., Lauer et al. 1992; Crane et 
al. 1993; Stiavelli, M0ller & Zeilinger 1993). This indepen- 
dently prompted Dehnen (1993, hereafter D93), Tremaine et 
al. (1994, hereafter T94) and Carollo (1993; hereafter C93), 
to study generalizations of the Jaffe and Hernquist models 
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with p oc r 7 at small radii. These '7-models' have a lumi- 
nous mass density of the form: 



Pl(t) 
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with a a scale-length, Ml the total luminous mass of the 
system and < 7 < 3. The profile falls off proportional 
to r -4 at large radii. The models of Jaffe and Hernquist 
are recovered for 7 = 2 and 7 = 1, respectively. The gravi- 
tational potential associated with the luminous density (1) 
is obtained from the Poisson equation, and has the simple 
form (C93; D93; T94): 
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where G is the gravitational constant. 

Many dynamical properties of the 7-models with an 
isotropic DF, / = f(E), were given by D93 and T94. Here 
we discuss the dynamical properties, and particularly the 
VPs, of 7-models with DFs of the anisotropic Osipkov- 
Merritt (079; M85) type. Our study is motivated by the 
fact that only few detailed VP studies for anisotropic sys- 
tems are available in the literature (e.g., Dejonghe 1987; 
Dehnen & Gerhard 1993; Gerhard 1993; van der Marel & 
Franx 1993). 

In Section 2 we show that the VPs of models with an 
Osipkov-Merritt DF can be written as a one-dimensional 
integral. This integral is no more complicated than the cor- 
responding integral for the isotropic case (included in the 
Osipkov-Merritt models as a limiting case). In Section 3 we 
discuss the range of anisotropy radii for which the Osipkov- 
Merritt 7-models are physical and likely to be stable, and 
we present the VPs for models with and without dark ha- 
los. To facilitate the comparison with observational data 
we characterize the predicted VP shapes by means of their 
Gauss-Hermite moments (e.g., van der Marel & Franx 1993). 
We study the variation of the Gauss-Hermite moments as a 
function of projected radius, cusp steepness 7 and (radial) 
anisotropy of the model. Concluding remarks follow in Sec- 
tion 4. 

Most of the results presented here are contained in C93. 
While work on the present manuscript was in progress, two 
papers appeared by Hiotelis (1994a, b). He studied some 
properties of Osipkov-Merritt DFs for models with the Jaffe 
luminosity density (7 = 2), both with and without a dark 
halo. He briefly discussed the VPs of the models, calculated 
by 'brute force' numerical evaluation of three-dimensional 
integrals. Our models are more general than his, in that we 
consider a range of values for 7 and discuss in more detail 
the dependence of the observables on the model parameters. 



2 OSIPKOV-MERRITT MODELS 

We derive a general expression for the VPs of Osipkov- 
Merritt models, valid for an arbitrary luminous mass density 
Pl{t), and for both self-consistent and non self-consistent 
models. In the former case only the luminous mass density 
contributes to the gravitational potential: $(?•) = Ql{t)- 
In the latter case there is an additional contribution from 



a dark matter component: $(?•) = Ql{t) + ^d{t)- Fol- 
lowing BT we define the relative potential of the system as 
^(t) = —$(?•). The relative potential at infinity, Vt<x>, can 
either have a finite value (such as in a system of finite mass) 
or be —00 (such as in a system embedded in a logarithmic 
potential). In the former case we set to zero by adding 
a suitable constant of integration, without loss of generality. 

We employ two coordinate systems for our models. Pro- 
jected dynamical quantities are described with cylindrical 
coordinates (R, <t>,z), with R the projected radius on the 
sky, <f> the azimuthal angle on the sky and z along the line of 
sight. Internal dynamical quantities are described with the 
corresponding spherical coordinates {r,d,<f>). 

2.1 Distribution functions 

We consider the special class of non-rotating models for 
which the DF is of the Osipkov-Merritt type (079; M85; 
BT): 

f(E,L 2 ) = f(Q), with Q = E - (L 2 /2r 2 a ), 



(3) 



and the additional constraint that 



f(Q) = for Q < (4) 

Here E = * - (v 2 /2) is the 'relative energy' per unit mass 
and L 2 = r 2 (y1 + v 2 ^) is the squared angular momentum 
per unit mass. In a system with finite total mass, E is the 
binding energy of a star per unit mass, and the constraint (4) 
implies that there are no unbound stars in the system and 
no stars on certain high angular momentum bound orbits. 
The velocity dispersion anisotropy as measured by Binney's 
(1980) /3-parameter is (M85): 



P(r) = 1 



(5) 



2{vl) rl + r 2 ' 

The free parameter r a in the definition of Q is therefore 
called the anisotropy radius. Inside r a the velocity distribu- 
tion is nearly isotropic, while outside r a it becomes increas- 
ingly more radially anisotropic. For fixed pl, the unique 
isotropic model with / = f(E) is recovered in the limit 
r a — > 00. 

The DF for the Osipkov-Merritt models follows from 
the luminous mass density and the potential through an 
Abel transform, similar to the one in Eddington's formula 
for f(E) in the isotropic case (BT): 

1 dG(Q) 



/((3) -2W2 dQ ' 
where the function G{Q) is denned as: 
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and we have written 
r 2 

pq(r) = (1 + -j)pl(t). 
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For the self-consistent 7-models with pl as in equation (1), 
G{Q) is elementary for special values of 7 (H90; D93; T94). 
In the general case, G{Q) is most easily computed using the 
radius rather than the potential as the independent variable, 
i.e.: 

dpQ dr 



G(Q) 
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(9) 
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The inversion r = ^(vt) required to obtain the lower bound of 
the integral must generally be evaluated numerically. It can 
be evaluated analytically for the self-consistent 7-models 
with potential (2): 



■/{[ 



(7 = 2), 
l}. (7/2). 



The DF f(Q) given by equation (6) can be evaluated for any 
luminous mass density pl{t), potential ^(r) and anisotropy 
radius r a . It is physical when the function G is a monoton- 
ically increasing function of Q. 

2.2 Velocity dispersions and total kinetic energy 

The velocity dispersions {v 2 } and (vg) = (■u 2 4 ) = \{v\_) as a 
function of radius in Osipkov-Merritt models can be com- 
puted by integration of v 2 f(Q) and Vgf(Q), respectively, 
over all velocities. Alternatively, they can be found by in- 
tegration of the single Jeans equation for spherical models 
with / = f(E, L 2 ), given in eq. (4-30) of BT. The result is 
(cf. M85): 
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Calculation of pl{v 2 } can be done conveniently by trans- 
forming to r = ^(vt). For the self-consistent 7-models the 
integral for pl{v 2 } is elementary when 47 is an integer. 

The total kinetic energies T T and Tj_ in the radial and 
tangential directions, respectively, are denned as 



00 

T r = 2-k J dr r 2 pl{vI), T± = 2t J dr r 2 Pl{v 2 ±). (12) 



The ratio 2T r /T± equals 1 in an isotropic model, and is 
larger than one in f(Q) models with finite r a . Its value is 
a global measure of the dominance of radial over tangential 
motion. Models with too much radial motion are likely to be 
unstable to triaxial deformation by means of the radial or- 
bit instability (Polyachenko 1981; Fridman & Polyachenko 
1984). A useful rule of thumb (although not an exact in- 
stability criterion) is that 2T r /T± should not exceed some 
(model dependent) threshold value (e.g., Barnes 1985; Mer- 
ritt & Aguilar 1985; Barnes, Goodman & Hut 1986; Merritt 
1987; Aguilar 1988). 

Substitution of expressions (11) in the definitions (12) 
and exchange of the order of integration leads to 
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where W is the total potential energy 



W = ^\ drr 1 



dr 
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It follows that T r +T± = —W/2, which is the virial theorem. 
For a self-consistent model d^(r)/dr = —GM(r)/r 2 , where 
M[r) = 47r dr r 2 Pl{t) is the (luminous) mass enclosed 
inside radius r. Both M(r) and W can be given explicitly 
for the 7-models (D93; T94): 



(1°) M(r) 
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W is infinite when 7 > 5/2. We study the behaviour of 
2T r /T± for the Osipkov-Merritt 7-models in Section 3.1. 

2.3 Evaluation of the velocity profiles 

By definition, the velocity profile VP is obtained by inte- 
grating the DF f(E, L 2 ) over the line of sight and over the 
two velocity components that lie in the plane of the sky. 
Hence, the normalized VP is 



VP(fl, v z ) = ^ j ' dz J J dv R dv* f(Q), 
where S is the projected surface density: 
E(R) = J dz JJJ dv R dv<t, dv z f(Q) 



Jr 



(17) 
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The VP does not depend on the azimuthal angle <f> on the 
plane of the sky, because of the spherical symmetry. We 
now show how the triple integral (17) can be simplified. 

Using (vr, v$, v z ) as coordinates in velocity space, the 
squared angular momentum is: 

L 2 = r 2 [(vr cos 6 — v z sin 6) 2 + v^] ■ 



(19) 



Hence, 



(20) 
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At a fixed point in space, Q is quadratic in the velocities, 
and can be cast in diagonal form by replacing (vr, v<p) by 
new velocity space coordinates (cj,£), defined through 



r" cos » sin » / r a + r 

VR = — ; — V z + 



r 2 + r 2 cos 2 1 
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v$ = ui cos £. 

Substitution in equation (20) yields: 



r 2 - R 2 
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which does not depend on £. We therefore transform the 
integral over dvR dv$ in equation (17) into an integral over 
d£ dui: 
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We carry out the trivial integration over £, and transform 
the integration over uidui to an integration over Q, by means 
of (22). With equation (4) this gives 

f f /"Qmax 

JJf(Q)dv R dvt = 2-Kg(r,R)J f(Q)dQ, (24) 



where 

g{r, R) ■■ 
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The upper limit Qmux follows from (22) upon substitution 
of u = 0: 



<?max(7-, R, v\) = 



+ r 2 



R? 2 



(26) 



Substitution of expression (24) in equation (17) yields for 
the VP: 

where, as in equation (18), the integral over dz is trans- 
formed to an integral over dr. With equations (6) and (7) 
this reduces to 
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where I(R, v z ) is the set 

I(R, v z ) = {r : r > R A Q^^{r, R, v 2 z ) > *<x>}. 



(29) 



In equation (28) we have used the fact that G(Stoo) = 
in models with finite total luminous mass. Once G(Q) has 
been calculated and tabulated as a function of Q using equa- 
tion (9), expression (28) allows the VPs to be calculated as 
one-dimensional integrals for all (R, v z ) combinations, with- 
out the explicit evaluation of the DF. The isotropic case is 
recovered in the limit r a —> oo, so that Q —> E, g(r, R) —> 1 
and Q max -> - v\j2. 

If = — oo, then I(R,v z ) is the interval [R, oo) and 
the VP extends to infinitely large velocities. In the follow- 
ing we restrict ourselves to models with finite total mass, 
and hence Vtoo = 0. At given R, the maximum absolute 
line-of-sight velocity value \v z \ lna . K {r, R) reachable by the 
stars at a given intrinsic radius r is set by the requirement 
<?max(' - , R, v 2 z ) = 0. Hence, 



u z max 



(r,R) 



2V{r 2 + r 2 — R 2 ) 
rl + r 2 



(30) 



The VP(i?, v z ) extends over a finite range of velocities, 
—Vt(R) < v z < Vt(R), where Vt(R) is the 'terminal velocity', 
which is obtained by maximizing | , u z | m ax(' - , R) as a function 
of r on the interval [R, oo). At large radii, r/R —> oo, one 
has |'U z | m ax(' - , R) ~ •^/2$(7 - ), which is a decreasing function 
of radius. As a result, the maximum is attained at a finite 
value of R, i.e., either at r = R, or at a radius r for which 
d|^|ma* (r, R)/dr = 0: 

-(rl + r 2 - J R 2 ) + ^4 =0 . (31) 



dr 



+ r z 



In an isotropic model r a —> oo, and equation (30) re- 
duces to |'U z | m ax(' - , R) = y^2$(7^). Over the interval [R, oo) 



this attains its maximum at r = R. In an isotropic model 
the terminal velocity is thus determined by those stars along 
the line of sight that are at the tangential point (deepest in 
the potential well), have zero energy and have their velocity 
vector along the line-of-sight. The angular momentum of 
these stars satisfies L 2 = 2i? 2 Vt(.R). In anisotropic models 
with finite r a there are no such stars, because of the imposed 
constraint (4). It follows that the terminal velocity need not 
necessarily be attained by stars at the tangential point. In 
addition, the terminal velocity is lower than in the isotropic 
model. 

For anisotropic models one must search (numerically) 
for solutions of equation (31). We found that for the 7- 
models there is generally either no, or only one solution. 
When no solutions are found, the terminal velocity is at- 
tained at r = R. When a solution ro > R is found, it is 
attained at the radius r = ro. At values of v z other than the 
terminal velocity, the set I(R,v z ) defined by equation (29) 
generally reduces to an interval [Vi,^], with R < n < r2 < 
00. The values of n and 7^2 must be calculated numerically 
for each combination of (R,v z ). 

The VPs to be discussed in the following sections were 
calculated by numerical integration of equation (28) for an 
array of v z values spaced linearly between —Vt(R) and Vt(R). 



3 RESULTS FOR OSIPKOV-MERRITT 
7-MODELS 

First we determine which combinations of values of 7 and r a 
correspond to physical distribution functions f(Q) > 0, and 
we discuss which of these models are expected to be dynam- 
ically stable. Then we investigate the kinematic properties 
of the Osipkov-Merritt 7-models. We consider the self- 
consistent case and also investigate the effect of adding the 
gravitational potential of a dark halo to the potential (2). 
We restrict ourselves to models with 7 < 2. 

3.1 Existence and stability 

The isotropic self-consistent 7-models (r a —> 00) all have 
physical distribution functions f(E) (e.g., T94). When r a 
decreases, the fraction of stars on radial orbits in an f(Q) 
model increases. In the limit r a —> only radial orbits are 
populated. For such a radial orbit model to be physical, the 
density must diverge as r~ 2 , or steeper, in the centre (BT). 
Thus, f(Q) for the anisotropic 7-models with 7 < 2 must 
become negative for sufficiently small values of r a in order 
to compensate for the overpopulation of the centre by the 
radial orbits needed to reproduce the density in the envelope 
of the model. We have numerically determined the area in 
the (7, 7 , a )-plane where the Osipkov-Merritt 7-models are 
physical. The thick drawn solid curve in Figure 1 marks the 
lower boundary of the physical self-consistent models. As 
expected, the minimum allowed value ofr a /a increases with 
decreasing 7, and ranges from at 7 = 2 to 0.44 at 7 = 0. 

N-body experiments have shown that f(Q) models 
with 7 = 2 are prone to the radial orbit instability when 
2T r /T± > 2.5 or, equivalently, r a < 0.25a (Merritt & 
Aguilar 1985). Little is known about the stability of 
Osipkov-Merritt 7-models with 7 < 2. Antonov's (1962) 
sufficient condition for stability against radial perturbations 
(Dejonghe & Merritt 1988) can be applied easily to the 
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Figure 1. The (7, r a )— plane for Osipkov-Merritt 7— models. 
Models in the hatched area below the thick solid curve have 
f(Q) < for some values of Q, and are unphysical. The dashed 
curve is r a = 0.25r^, with th the half— mass radius. This curve 
roughly represents the lower boundary of the region that corre- 
sponds to models which are not affected by the radial orbit in- 
stability. The thin solid curves are contours of constant 2T r /T±. 
The contour values are indicated. 

Osipkov-Merritt 7-models, but we find that it is inconclu- 
sive: it confirms that the isotropic models (r„ —> 00) are 
stable, but gives no such guarantee for finite r a . 

The half-mass radius th of the 7-models is given by 
Th /a = l/(2 1/(3_7) - 1), so that r h = a when 7 = 2. 
To good accuracy i? e ff = 0.75th, where i? e ff is the radius 
which contains half the projected surface density of the 7- 
models (D93). A naive generalization of the stability con- 
dition r a > 0.25a when 7 = 2 is to take r a = 0.25rh as 
the stability boundary. This relation is shown as the dashed 
line in Figure 1. However, the relative importance of radial 
orbits in models with fixed r a /rh but different 7 increases 
when 7 decreases. This causes the steeper slope of the con- 
tours of constant 2T r /T± in Figure 1 (thin solid lines). If the 
stability boundary is given by 2T r /T± = 2.5 also for f(Q) 
models with 7 < 2, then models with r a /<i < 1.6 — 0.677 are 
radial orbit unstable. 

The determination of the precise location of the stability 
boundary in the (7, r^-plane requires a normal mode anal- 
ysis (e.g., Weinberg 1991), or careful N-body experiments. 
In what follows, we will only discuss kinematic properties 
of physical Osipkov-Merritt 7-models which are likely to be 
stable. 

3.2 VPs for the self— consistent models 

Since our models have zero mean streaming, their VPs are 
symmetric. In order to facilitate the comparison with obser- 
vational data, we express the VP as an orthogonal Gauss- 
Hermite series (e.g., van der Marel & Franx 1993; Gerhard 
1993) with parameters: (i) the line strength and the disper- 
sion a of the best-fitting Gaussian VP; and (ii) the dimen- 



sionless Gauss-Hermite moments zi,z$, . . . that measure de- 
viations from the best-fitting Gaussian with zero mean (van 
der Marel et al. 1994b). A value zi > generally indicates 
that the VP is more centrally peaked than a Gaussian, a 
value zi < that it is more flat-topped. To study the ob- 
servable properties of our models, we consider only a and zi 
since the signal-to-noise ratio and instrumental resolution 
of spectroscopic data are often insufficient to measure Z6 and 
higher order coefficients. Typical uncertainties in published 
Zi measurements for elliptical galaxies are of the order of 
0.02. 

Figure 2 illustrates the observable kinematic properties 
of self-consistent Osipkov-Merritt models with 7 = 1. In 
this case the half-mass radius th and the scale length a 
are related through th = 2.4142a, and the effective radius 
Reg = 1.815a (D93). The panels on the left show the ef- 
fect of anisotropy on the radial a and Zi profiles, while the 
panels on the right show the normalized VPs at R = a 
and R = r a , as a function of v z /cr. Plotted are the mod- 
els with r a = 0.3th (dotted), r a = 0.5th (short-dashed) 
T a = th (long-dashed) and r a = 2th (dot-dashed). The 
solid curve has r a —> 00, and corresponds to the isotropic 
model. The model with r a = 0.3th lies above the dashed 
line in Figure 1, but below the contour 2T r /T± = 2.5, and 
its stability is not guaranteed. All VPs (including those 
for the isotropic model) differ appreciably from a Gaussian, 
and show a strong variation with radius. The a profiles of 
the anisotropic models lie below the isotropic curve at large 
radii, while they lie above it at small radii. When R 3> r a 
the motion is strongly radial, and only a small fraction of 
this is in the direction of the line of sight. When R < r a , cr 
is large because of the contribution of fore- and background 
stars at large galactocentric distances, which now have a 
substantial part of their motion directed towards the ob- 
server. As expected, the anisotropic models have ^-profiles 
that are similar to the isotropic one when R <C Ta- When 
R ~ r a the Z4-values are negative, because stars at radii r 
slightly larger than r a , though less in number, have a larger 
part of their motion along the line of sight, and hence con- 
tribute strongly to the VP at large velocities. This effect 
is stronger when r a becomes smaller than a, because then 
the density profile in this radial range falls off like r ', and 
the typical local rms velocities are still increasing outwards 
(■u cc r°' b when p cc 7 --1 ). When R 3> r a the anisotropic z$- 
profiles increasingly deviate from the isotropic curve towards 
more positive values at larger radii due to the increasing ra- 
dial anisotropy. Now most of the motion is perpendicular 
to the line of sight, and hence most of the contribution to 
the VP is at small velocities, so that the VP is much more 
peaked than a Gaussian. 

Similar results are obtained for other values of 7. In 
Figure 3 we illustrate the case 7 = 3/2 (for which th = 
1.7024a and i? e ff = 1.276a; D93). Again the stability of the 
model with r a = 0.3th is not guaranteed. The cusp in the 
density profile now has a steeper slope, so the stars in the 
centre contribute more strongly to the projected kinematic 
quantities. 

Figure 4 shows the effects of varying the slope 7 of the 
inner density profile for the anisotropic self-consistent model 
with r a = a. Plotted are the velocity profiles at R = a and 
R = 0.2a, and the radial a and zi profiles for 7 = 0.5 (long- 
dashed; R e f[ = 2.358a), 7 = 1.0 (short-dashed), 7 = 1.5 
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Figure 2. Projected kinematic profiles for the self— consistent 7 = 1 (H90) model described in the text. Different curves correspond to 
different values for the ratio r a /a of the anisotropy radius and the length scale of the mass density: r a = 0.72a (dotted curve), r a = 1.21a 
(short— dashed curve), r a = 2 Ala (long— dashed curve), r a = 4.83a (dot— dashed curve) and r a — > oo (solid curve; isotropic model). These 
curves correspond to values of the ratio r a /r^ = 0.3,0.5, 1.0, 2.0, and oo, where is the half— mass radius of the mass density, a) radial 
profiles of the projected velocity dispersion a in units of <\J ~GM~j J /a. b) radial profiles of the fourth Gauss-Hermite coefficient Z4 . A value 
Z4 > generally indicates that the VP is more centrally peaked than a Gaussian, a value Z4 < that it is more flat-topped, c) normalized 
line-of-sight velocity profiles as function of v z /a, at projected galacto centric distance R = a. d) idem, at R = r a - The velocity profiles 
are symmetric, and only the positive velocity part is shown. 




Figure 3. Projected kinematic profiles similar to Figure 2, but now for the self— consistent 7 = 3/2 model. As in Figure 2, the curves 
are for the values 0.3, 0.5, 1.0, 2.0, and 00 of the ratio r a /r^. These now correspond to r a = 0.51a, 0.85a, 1.70a, 3.40a and 00. 
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Figure 4. Projected kinematic profiles for self— consistent 7— models with anisotropy radius r a = a. The different curves correspond to 
different values for the cusp steepness 7: 7 = 0.5 (long— dashed curve), 7 = 1.0 (short— dashed curve), 7 = 1.5 (dotted curve) and 7 = 2.0 
(solid curve). Panels a) and b) show the a and Z4 profiles, similar to Figure 2. Panels c) and d) show the velocity profiles at projected 
galactocentric distances R = a and R = 0.2a, respectively. 



(dotted) and 7 = 2.0 (solid; R eff = 0.744a). The model 
with 7 = and r a = a is not shown since it is possibly 
unstable (cf. Figure 1). Varying 7 does not affect the dy- 
namical structure at large radii, so the profiles coincide for 
R ^ T a . When 7 increases towards 2, i.e., for steeper central 
density profiles, the projected properties of the models are 
increasingly dominated by the stars in the central region. 
Here the velocity distribution of the f(Q) models is nearly 
isotropic. For 7 — > 2 (as in the singular isothermal sphere), 
it approaches a Gaussian. As a result, increasing 7 to 2 de- 
creases the Zi values inside R = r a towards 0. As expected, 
the asymptotic behaviour of a for R —> is similar to that 
in the isotropic models (see also D93 and T94). 

3.3 VPs for models with dark halos 

We now consider the same models embedded in a dark halo. 
We assume that the dark potential is also represented by a 
7-model, but with different scale-length 3d > 3 and total 
mass Mb than for the luminous matter. We only investigate 
the case where the value of 7 is the same for the luminous 
and the dark mass density. Our choice of dark halo potential 
results in a nearly flat rotation curve in the observationally 
accessible range when an /a is sufficiently large, and has the 
advantage of giving a finite total mass. 

The isotropic models with dark halos are physical for 
all values of 7. We have found that the location of the 
lower boundary in the (7, 7 , a )-plane of the physical 7-models 
embedded in a dark halo with a range of masses Mb I Ml 
and scale-lengths an I a differs insignificantly from the thick 
solid curve in Figure 1. We expect that the location of the 
stability boundary for the models with dark halos also does 



not differ very much from the self-consistent boundary (cf. 
Stiavelli & Sparke 1991). 

The VPs were again calculated as described in Sec- 
tion 2.2. The only differences from the self-consistent case 
are that now = — $1 — $d, and that the inversion ^(vt) 
required for equation (9) must be evaluated numerically. 
Since the total mass of the halo is finite, we can again take 
VP 00 = 0, and G(ty oo) = 0. This means that also for these 
non self-consistent models the calculation of the VP is done 
without an explicit evaluation of the DF. 

Figure 5 illustrates the effect of adding a dark halo on 
the behaviour of the VPs at R = a and R = 0.2a, and 
of a and zi as a function of radius for the isotropic mod- 
els with 7 = 3/2. Plotted are the models with an = 5a 
and M D /M L = 10.0 (short-dashed), Mb /Ml =5.0 (dot- 
ted) and Mb I Ml = 0.0 (solid). The major effect of the 
introduction of the dark halo is the added depth of the po- 
tential well, which results in increased velocity dispersion at 
all radii. The velocity distribution remains isotropic, and 
hence the zi values hardly change. 

Figure 6 shows the similar results for the anisotropic 
models with r a = a, again for 7 = 3/2. Again the a values 
increase, as expected. The Zi profiles of the models with 
dark halos lie systematically above the no dark halo curve, 
but the differences are not large. The dip in Zi seen in 
the models without dark halo (Figure 3) is less pronounced. 
This is caused by the fact that the extended dark halo po- 
tential softens the radial variation of the typical orbital ve- 
locities so that the contribution to the VP of the stars at 
radii r > r a is diminished. 

Figure 7 shows the results for the same anisotropic mod- 
els, but now with Mb = 10 Mi,, for different ratios of the 
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Figure 5. Projected kinematic profiles for isotropic models (r a — > oo) with the luminous mass density of a 7 = 3/2 model, with total 
luminous mass Ml and scale-length a, embedded in a dark halo that also has the mass density of a 7 = 3/2 model, but has total dark 
mass M.j) and scale-length . The different curves correspond to different values for the ratio /Ml ■ Mp /Ml = 10.0 (short— dashed 
curve), Mj) /Ml = 5.0 (dotted curve) and Mp /Ml = 0.0 (solid curve). All models have = 5a. The panels a) and b) show the a (in 
units of <\J GMl /cl) and Z4 profiles, similar to Figure 2. Panels c) and d) show the velocity profiles at projected galacto centric distances 
R = a and R = 0.2a, respectively. 




Figure 6. Projected kinematic profiles for 7 = 3/2 models with r a = ct, embedded in a 7 = 3/2 dark halo. The scale— length of the 
dark halo is taken equal to 5a in all models, while the ratios Mp /Ml are varied as in Figure 5. 



dark to luminous scale-lengths ac/a = 10,5,2. The self- 
consistent model (Mc = 0) is shown for comparison (solid 
curve). The behaviour of the z± profiles is similar as in Fig- 



ure 6, and this has the same reason: when ac/a increases, 
the radial variation of the potential is increasingly softened, 
and the dip in the z± profile is smoothed out. The a profiles 



Velocity profiles of Osipkov-Merritt models 9 





Figure 7. Projected kinematic profiles for 7 = 3/2 models with r a = a, embedded in a 7 = 3/2 dark halo. This time all models 
have /Ml fixed at the value 10, and the ratio ajj /a is varied. The different curves correspond to ajj /a = 10.0 (long— dashed curve), 
a D l a — 5.0 (short— dashed curve), ajj /a = 2.0 (dotted curve). The model with Mjj = is shown for comparison (solid curve). 



with dark halo are all for the same total mass, and hence 
they converge at radii R 3> an ■ The differences at smaller 
radii are caused by the increasing central concentration of 
the potential when an decreases. 

4 CONCLUDING REMARKS 

The Osipkov-Merritt spherical dynamical models have the 
special property that their distribution functions are a func- 
tion of one integral of motion: / = /(<?), where Q is a 
linear combination of the binding energy E and the square 
of the total angular momentum, L 2 . As such, they are the 
simplest spherical models with the attractive property that 
the intrinsic velocity distribution is isotropic in the central 
region, and becomes increasingly radially anisotropic going 
outwards. This behaviour is in good qualitative agreement 
with dissipationless collapse models for the formation of el- 
liptical galaxies (e.g., van Albada 1982; Bertin & Stiavelli 
1993). 

We have shown that the 7-models introduced by D93, 
C93 and T94 correspond to physical distribution functions 
f(Q) > in a large region of the (7, r^-parameter space: 
only the hatched area in Figure 1 is excluded. All models 
with 7 and r a above the dashed line in Figure 1 are expected 
to be dynamically stable, but this conclusion needs to be 
confirmed by detailed stability analysis. 

The calculation of the line-of-sight projected veloc- 
ity distribution (the VP), which generally requires a triple 
numerical integration, reduces to a single quadrature in 
/ = f{Q) models, and therefore is no more difficult than for 
isotropic / = f(E) models. By contrast with the isotropic 
models, however, the VPs of Osipkov-Merritt models gen- 
erally show a strong variation with projected galactocentric 



radius. We have illustrated this range of VP behaviour for 
the 7-models. The VPs may be more peaked than a Gaus- 
sian at small radii, significantly less peaked at intermediate 
radii, and again more peaked at large radii. This behaviour 
is caused by the trade-off between the increasing radial bias 
at large intrinsic radii, and a decreasing contribution of the 
radial component of the velocity ellipsoid to the line-of-sight 
velocity at larger projected radii. We quantified the VP 
variations by the dimensionless shape-parameter Zi, which 
can vary by as much as 0.1-0.2. This is an order of mag- 
nitude larger than the typical measurement uncertainties in 
Zi. These VP variations are thus detectable with current 
observational techniques. 

The spherical 7-models provide a good fit to the surface 
brightness profiles of round elliptical galaxies. We expect 
that the VP properties described here for 7-models with 
and without dark halos will be useful for the interpretation 
of high quality kinematic data for such galaxies. 
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